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Abstract 

Direct  electrochemical  oxidation  (DECO)  solid-oxide  fuel  cells  (SOFCs)  offer  the  potential  to  generate  electrical  power  from  hydrocarbon 
fuels  without  the  need  for  upstream  fuel  processing,  such  as  reforming.  However,  with  pure  hydrocarbon  fuel  entering  the  flow  channels 
at  temperatures  typically  above  700  °C,  fuel  pyrolysis  can  cause  molecular- weight  growth  and  the  formation  of  deleterious  carbonaceous 
deposits.  This  paper,  which  develops  a  plug-flow  model  for  fuel  (natural  gas  surrogate)  within  the  anode  channels,  considers  the  elementary 
gas-phase  chemical  kinetics  of  fuel  pyrolysis  and  oxidation.  It  also  considers  the  limiting  case  of  local  chemical  equilibrium.  Formation  of 
cyclic  hydrocarbon  species  is  used  to  indicate  deposit  propensity.  Results  show  that  the  likelihood  of  deposit  formation  depends  strongly 
on  cell  temperature,  current  density,  and  residence  time.  Generally  speaking,  equilibrium  favors  deposit  formation  early  in  the  channel 
whereas,  owing  to  limited  residence  time,  the  homogeneous  finite-rate  kinetics  predicts  relatively  low  levels  of  deposit  precursors.  In  the 
downstream  portions,  because  of  electrochemical  oxygen  flux  though  the  electrode-electrolyte  membrane,  chemical  equilibrium  shifts 
strongly  away  from  deposit  formation  to  volatile  carbon-oxygen  species.  However,  the  homogeneous  finite-rate  kinetics  predictions  show 
a  continuing  increase  in  coking  propensity. 
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1.  Introduction 

Depending  on  the  particulars  of  electrode  and  elec¬ 
trolyte  materials,  direct  electrochemical  oxidation  (DECO) 
solid-oxide  fuel  cell  (SOFC)  systems  typically  operate  at 
temperatures  between  700  and  900  °C  [1].  With  a  pure 
hydrocarbon  gas  entering  the  anode  flow  channels  at  these 
temperatures,  it  is  well  known  that  the  fuel  can  pyrolyze. 
These  conditions  can  lead  to  molecular-weight  growth  and 
ultimately  to  the  formation  of  cyclic  and  polycyclic  aromatic 
hydrocarbons,  which  are  indicators  of  a  propensity  to  form 
carbonaceous  deposits.  These  deposits  can  damage  porous 
electrode  structures,  thereby  diminishing  cell  life  and  per¬ 
formance.  Deposit  formation  depends  on  temperature  and 
residence  time  as  well  as  the  local  chemical  composition 
of  the  stream.  Pure  fuel  enters  the  anode  channels  and 
is  electrochemically  oxidized,  producing  E^O  and  CO2. 
Increasing  product  concentrations  toward  the  downstream 
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end  of  the  anode  channel  possibly  inhibit  fuel  pyrolysis  by 
altering  residence  time  and  reaction  rates. 

The  principal  objective  of  this  work  is  to  character¬ 
ize  the  gas-phase  chemistry  along  the  anode  channel  of  a 
DECO-SOFC  as  a  function  of  cell  design  and  operating 
conditions.  The  cell  current  density  determines  the  rate 
at  which  hydrocarbon  fuel  is  electrochemically  converted 
to  H2O  and  CO2.  Along  with  current  density,  the  fuel 
flow  rate  and  channel  dimensions  affect  the  residence  time 
and  chemical  composition,  thus  altering  the  likelihood  of 
molecular-weight-growth  chemistry. 

The  homogeneous  kinetic  model  considers  a  large,  el¬ 
ementary  reaction  mechanism  that  is  well  suited  to  un¬ 
derstanding  fuel  stability  and  molecular-weight  growth  for 
gas-phase  hydrocarbon  systems.  In  particular,  we  look  for 
the  production  of  cyclic  hydrocarbon  molecules  as  indica¬ 
tors  of  deposit  formation,  much  as  these  species  are  used  as 
indicators  of  soot  production  in  combustion  processes  [2]. 
Equilibrium  calculations  are  performed  along  the  channel, 
providing  a  basis  for  comparison  between  the  kinetic  pre¬ 
diction  and  local  thermodynamic  equilibrium. 
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In  addition  to  gas-phase  reactions,  heterogeneous  kinetics 
play  a  significant  role  in  electrochemical,  oxidation  and  py¬ 
rolysis  reaction  pathways.  However,  due  to  the  complexity 
of  DECO-SOFC  systems,  we  have  focused  on  the  homoge¬ 
nous  kinetics  in  this  work.  This  is  done  in  part  to  identify 
the  relative  magnitude  of  homogeneous  fuel  chemistry  and 
the  effects  of  electrochemical  oxidation.  Local  equilibrium 
predictions  provide  some  indication  of  what  might  happen 
to  the  fuel  stream  in  the  presence  of  catalytically  accelerated 
reactions. 


1.1.  Kinetic  and  equilibrium  models 

Koh  et  al.  [3]  presented  an  equilibrium  model  for  solid 
carbon  deposition  in  SOFC  systems  using  dry  hydrocarbon 
fuels.  The  relationship  between  current  density  and  the  flux 
of  oxygen  atoms  into  the  system,  along  with  the  design  and 
operating  conditions  of  the  cell,  formed  the  basis  for  a  map  of 
equilibrium  solid  carbon  as  a  function  of  system  parameters. 
However,  this  work  did  not  account  for  finite-rate  kinetics  in 
the  fuel  stream.  Consideration  of  both  kinetics  and  chemical 
equilibrium  is  essential  to  gain  an  improved  understanding 
of  the  fuel  chemistry  in  DECO-SOFC  systems.  The  kinetic 
model,  when  combined  with  a  fluid-flow  model,  predicts 
spatial  variation  in  gas-phase  species  composition  along  the 
channel.  A  free-energy  minimization,  given  the  temperature, 
pressure  and  local  atomic  species  distribution  from  the  so¬ 
lution  to  the  kinetic  model,  gives  the  chemical  equilibrium 
state  of  the  local  mixture  at  a  particular  point  in  the  channel. 


1.2.  DECO-SOFC  background 


DECO  type  SOFC  systems  are  a  rapidly  advancing  tech¬ 
nology,  potentially  offering  a  highly  efficient  alternative  for 
electrical  power  generation  from  hydrocarbon  fuels  with¬ 
out  upstream  fuel  preparation  such  as  reforming  or  water 
addition  [4].  Fig.  1  illustrates  some  of  the  essential  ele¬ 
ments  typical  of  a  planar  SOFC  architecture.  The  heart  of 
the  cell  is  the  electrochemically  active  membrane-electrode 
assembly  (ME A),  which  consists  of  a  dense  ceramic  elec- 


Fig.  1.  Illustration  of  a  segment  of  two  fuel-cell  layers,  including  the 
electrode-electrolyte  assembly  (MEA),  the  interconnect,  and  the  flow 
channels  for  fuel  and  air.  A  counter-flow  situation  is  illustrated  here,  but 
co-flow  and  cross-flow  configurations  are  also  common. 


trolyte  sandwiched  between  porous,  ceramic-metal  compos¬ 
ite  electrodes.  The  electrolyte  is  impervious  to  gas  flow,  but 
is  an  oxygen-ion  (O2-)  conductor.  Reactants  and  products 
are  transported  through  typically  rectangular  flow  channels, 
with  a  cross-sectional  area  in  the  order  of  0. 5-1.0  mm2.  The 
length  and  number  of  channels  varies,  depending  on  the  ac¬ 
tive  MEA  area  and  hence  the  power  level  of  the  stack  [5]. 
Usually,  several  tens  of  cell  layers  are  stacked  and  electri¬ 
cally  connected  in  series.  Because  an  individual  cell  layer 
produces  an  output  voltage  approximately  0.5  V  at  peak 
power  density,  the  stack  configuration  increases  the  operat¬ 
ing  voltage  to  more  practical  levels. 

In  a  methane-fed  DECO-SOFC,  there  are  two  principal 
electrochemical  reactions.  On  the  cathode  side  of  the  elec¬ 
trolyte,  O2  is  reduced, 

02  +  4e-^202-,  (Rl) 

producing  negatively  charged  ions  that  are  transported 
through  the  electrolyte  to  the  anode  interface.  At  the 
electrolyte-anode  interface,  electrochemical  oxidation  reac¬ 
tions  occur,  which  can  be  expressed  globally  for  methane  as 

CH4  +  402-  ^  2H20  +  C02  +  8e“,  (R2) 

depleting  fuel  and  producing  CO2  and  H2O.  While  homo¬ 
geneous  reactions  produce  other  species  that  are  known  to 
participate  in  electrochemical  oxidation  reactions  (e.g.  H2 
and  CO)  under  these  conditions  very  little  CH4  is  converted 
to  such  species  in  the  gas  phase.  Electrochemical  oxidation 
also  frees  electrons  from  the  charged  species,  producing  an 
electrical  current  through  an  external  circuit.  While  the  ac¬ 
tual  process  involves  many  reaction  steps  and  complex  het¬ 
erogeneous  electrochemistry  at  the  three-phase  boundaries, 
the  global  representation  is  sufficient  to  account  for  the  elec¬ 
trochemically  produced  product  species  on  the  gas-phase 
chemistry  in  the  flow  channels. 

Based  on  the  global  anodic  electrochemical  reaction  (R2), 
the  molar  fluxes  between  the  fuel  stream  and  the  MEA  of 
the  methane-SOFC  system  at  a  local  current  density  of  /e 
can  be  represented  as, 

Jcn4  =  -jp,  Jco2  =  -pp,  Ju2o  =  ~fp,  (1) 

where  F  is  the  Faraday  constant.  As  the  electric  cell  po¬ 
tential  Eceii  across  the  MEA  structure  can  be  assumed  to 
be  constant  along  the  channel  due  to  the  excellent  electric 
conductivity  of  the  interconnect,  the  current  density  /e  may 
vary  along  the  channel  as  a  result  of  fuel  and  oxidizer  con¬ 
sumption.  However,  the  relationship  between  the  current 
density  /e  and  the  electric  cell  potential  Eceii  is  complex 
[6].  It  depends  on  all  the  polarization  losses  (i.e.  concentra¬ 
tion,  activation,  and  ohmic  polarizations)  across  the  MEA 
structure,  thereby  depending  on  the  material  properties  and 
micro-structures,  and  operating  conditions,  particularly  on 
the  fuel  and  oxidizer  compositions.  To  incorporate  the  ef¬ 
fect  of  fuel  consumption  along  the  channel  into  the  current 
model,  the  following  simple  model  for  the  current  density  /e 
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is  used, 

U  =  ie[CH4F.  (2) 

where  i®  is  a  constant,  [CH4]  is  the  local  molar  concentration 
of  CH4  in  the  fuel  stream.  The  exponent  y  is  a  constant  de¬ 
termining  the  order  dependence  of  /e  on  CH4  concentration. 


2.  Model  structure 

The  Reynolds  number  (based  on  a  characteristic  channel 
diameter)  for  typical  anode-channel  flow  is  of  order  unity. 
Since  the  entry  length  for  small  Reynolds  numbers  is  short 
compared  to  the  channel  length,  we  make  a  plug-flow  ap¬ 
proximation  that  considers  only  the  mean  axial  velocity. 
Variations  in  temperature  and  composition  across  the  chan¬ 
nel  are  assumed  to  be  negligible.  This  assumption  was  jus¬ 
tified  by  a  boundary-layer  model  that  showed  negligibly 
small  radial  variation  in  composition  or  temperature  across 
the  channel  width.  We  chose  to  neglect  diffusive  transport 
in  the  axial  direction,  which  is  consistent  with  the  model 
developed  by  Iwata  et  al.  [7].  Furthermore,  only  isothermal 
flows  are  discussed  in  this  analysis,  although  temperature 
variation  would  affect  the  reaction  chemistry  and  impact 
cell  performance.  We  assume  nominal  atmospheric  pressure 
throughout  the  stream.  The  plug-flow  model,  which  is  based 
on  the  PLUG  software  [8]  that  is  in  the  Chemkin-III  pack¬ 
age,  has  been  extended  to  accommodate  the  electrochemi- 
cally  driven  flux  of  reactants  (CH4)  and  products  (H2O  and 
CO2)  that  enter  or  leave  the  channel.  Detailed  chemical  ki¬ 
netics  are  handled  through  the  Chemkin-III  [9]  and  Surface 
Chemkin-III  [10]  software  interfaces. 

2.7.  Governing  equations 


Fig.  2,  which  depicts  a  control  volume  inside  an  anode 
channel,  serves  as  the  framework  within  which  the  needed 
conservation  equations  are  derived.  The  control  volume  has 
differential  length  d z  along  the  flow  direction,  but  spans  the 
channel’s  width  and  height. 

Restricting  attention  to  steady,  one-dimensional,  plug 
flow,  the  mass  conservation  equation  is  expressed  as, 


A 


c 


d  (pu) 
d  z 


K  K 

=  «e  V  -4  Wk  +  «'  SkWk. 

k=  1  k=  1 


Interconnect 


Electrode¬ 

electrolyte 

assembly 


Electrochemical 

flux 


Fig.  2.  Illustration  of  a  control  volume  inside  the  fuel  channel  with  the 
H2O  and  CO2  flux  coming  through  the  anode  into  the  fuel  stream. 


The  species  molecular  weights  are  given  as  Wk  and  Sk  are 
the  molar  production  rate  of  gas-phase  species  by  surface 
reaction.  Although  the  results  presented  in  this  paper  neglect 
catalytic  chemistry,  the  terms  are  retained  so  that  the  model 
can  be  extended  to  include  surface  chemistry. 

In  Eq.  (3),  Ac  represents  the  cross-sectional  area  of  the 
channel.  The  parameter  a'Q  is  the  electrochemically  active 
area  per  unit  length  of  channel  (i.e.  the  perimeter  available 
for  electrochemistry).  In  Fig.  2,  a'e  would  be  the  width  of  the 
channel  floor  plus  the  width  of  the  active  surface  under  the 
interconnect  ribs.  The  parameter  a's  is  area  per  unit  length 
(i.e.  the  perimeter)  of  the  channel  wall  that  is  available  for 
heterogeneous  surface  chemistry. 

The  surface-chemistry  term  hWk  can  only  provide  a 
net  mass  loss  or  gain  to  the  channel  flow  if  there  is  deposi¬ 
tion  on  the  wall.  Alone,  catalysis  can  only  alter  the  chemical 
composition  of  the  flow,  but  cannot  alter  the  mass  flow  rate. 
It  should  be  noted  that  the  surface-chemistry  term  in  the 
overall  continuity  equation  vanishes  for  purely  catalytic  re¬ 
action  where  there  is  no  net  mass  exchange  with  the  surface. 
However,  the  individual  contributions  to  the  summation  are 
important  in  the  species  continuity  equations. 

For  each  species  k,  the  general  conservation  equation  may 
be  stated  as 

K 

JkWk  -YkJ2  JkWk 

k=  1 

+  a's  /  j* Wk-YkJ2skWk\+Acd>kWk.  (4) 

where  Yk  represent  the  mass  fractions,  cok  are  the  molar 
production  rates  by  homogeneous  reactions,  and  Wk  are  the 
molecular  weights. 

The  conservation  of  momentum  is  expressed  as 


The  wall  shear  stress,  rw,  is  represented  here  in  terms  of  the 
friction  factor  /  as 


Tw 

(1/2  W 


For  the  low-Reynolds  number  channel  flow,  the  product  Re  f 
is  a  constant.  For  rectangular  channels,  the  value  of  Re  f 
(based  on  the  hydraulic  diameter,  Df)  depends  on  the  chan¬ 
nel  aspect  ratio  a  =  77/  W  as 

(  3.4  \  puDh 

Ref  «  13.74+  10.38  exp - ,  Re  =  - — -.  (7) 

\  a  J  /a 


The  principal  role  of  the  momentum  equation  in  the 
plug-flow  analysis  is  to  determine  the  pressure  profile  along 
the  length  of  the  channel.  Assuming  that  the  thermodynamic 
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pressure  is  nearly  uniform  for  the  purposes  of  evaluating 
density,  chemical  and  thermodynamic  properties,  the  con¬ 
tinuity  equations  can  be  decoupled  from  the  momentum 
equation.  In  other  words,  for  a  given  pressure,  the  velocity 
and  species  profiles  can  be  determined  from  the  continuity 
equations  alone. 

An  ideal-gas  equation  of  state  is  used  to  evaluate  the  mass 
density,  p ,  in  terms  of  pressure  and  species  mass  fractions, 


RTYf=i  Yk/Wk' 

Because  the  flow  is  presumed  to  be  isothermal,  an  energy 
equation  is  not  required.  Nevertheless,  it  should  be  noted 
that  it  is  not  uncommon  in  practice  to  have  temperature 
variations  along  the  length  of  a  flow  channel.  The  details  of 
such  temperature  variations  depend  strongly  on  system-level 
considerations  like  the  stack  size  and  insulation  strategies. 
Since  the  present  model  is  concerned  with  a  single  channel,  it 
is  not  practical  to  incorporate  these  system-level  influences. 

2.2.  Elementary  gas-phase  mechanism 

The  detailed  chemical  mechanism  used  in  this  analysis 
is  an  extension  of  one  published  earlier  [11]  to  describe  the 
kinetics  of  C4  and  smaller  fuel  species.  It  has  been  extended 
to  selected  species  (n -hexane,  2, 3 -dimethyl-butane,  and 
cyclohexane)  by  allowing  for  both  dissociation  and  hydro¬ 
gen  abstraction  from  the  parent  fuel  species.  The  radicals 
formed  from  the  parent  in  these  processes  are  then  allowed 
to  undergo  beta-scission,  isomerization,  and  bimolecular 
reactions.  This  sequence  is  continued  until  the  resulting 
products  are  already  included  in  the  base  C4  mechanism. 
The  underlying  C4  mechanism  has  also  been  modified  to 
reflect  improved  understanding.  One  such  modification  re¬ 
flects  a  detailed  analysis  of  propargyl  radical  recombination. 
This  mechanism  is  substantially  similar  to  one  previously 
used  to  model  the  effect  of  temperature  on  soot  formation 
in  diesel  engines  [2]. 

The  mechanism  includes  many  pressure-dependent  re¬ 
action  rate  coefficients  and  the  estimates  for  atmospheric 
pressure  are  used  for  the  current  calculations.  Pressure 
dependencies  arise  in  those  reactions  in  which  a  new  chem¬ 
ical  bond  is  formed,  depositing  the  bond  energy  in  the 
newly-formed  molecule.  The  subsequent  fate  of  this  species 
depends  upon  the  relative  rates  of  the  various  uni-molecular 
reaction  pathways  (dependent  on  temperature)  to  the  col- 
lisional  stabilization  rate  (dependent  on  pressure).  The 
procedure  for  estimation  of  these  pressure-dependent  rate 
coefficients  for  the  various  product  and  stabilization  chan¬ 
nels  has  been  previously  described  [12]. 

Hydrogen-abstraction  reactions  included  in  the  mecha¬ 
nism  were  assigned  rate  coefficients  using  a  procedure  also 
described  previously  [13].  These  rate  coefficients  are  esti¬ 
mated  by  relating  the  parameters  for  a  specific  reaction  type 
(e.g.  abstractions  by  OH)  to  that  for  a  well  characterized 
hydrocarbon  analog  (e.g.  OH  +  C2H6).  Consistent  use  of 


such  reference  reactions  avoids  propagating  whatever  incon¬ 
sistencies  are  present  in  published  compilations.  Moreover, 
this  approach  accounts  for  the  different  temperature  depen¬ 
dencies  found  for  various  types  of  hydrogen  abstractions.  It 
is  now  well  recognized  [14]  that  many  hydrogen-abstraction 
reactions  show  upward  curvature  in  their  Arrhenius  plots, 
leading  to  appreciably  higher  rate  coefficients  at  high  tem¬ 
peratures  than  one  would  estimate  from  a  linear  extrapo¬ 
lation  of  low-temperature  measurements.  Such  temperature 
dependencies  have  to  be  included  in  order  to  predict  reli¬ 
able  rate  coefficients  under  high- temperature  conditions.  A 
convenient  way  to  include  this  dependence  is  to  express  the 
rate  constant  in  modified  Arrhenius  form, 


k  ~  A7^exp 


which  is  consistent  with  Transition  State  Theory. 

The  reaction  mechanism  proposed  by  Frenklach  and 
Warnatz  [15]  was  employed  for  describing  formation  of 
the  second  aromatic  ring  and  further  cyclizations  and  ox¬ 
idation.  The  ring-growth  sequence  is  continued  up  to  cy- 
clopentacoronene  in  this  model.  The  aromatic  rings  grow 
by  H-abstraction-C2H2-addition  reaction  sequence.  These 
reactions  were  included  to  provide  a  pathway  to  large 
polynuclear  aromatics  (PAH),  which  were  considered  to  be 
the  species  that  could  potentially  form  deposits. 

All  reactions  are  considered  reversible,  with  the  reverse 
rate  coefficients  computed  using  calculated  temperature- 
dependent  equilibrium  constants.  The  thermodynamic  data 
required  to  determine  the  equilibrium  constants  were  ob¬ 
tained  using  group  additivity  [16,17]  for  those  species  where 
experimental  thermodynamic  data  were  not  available. 

This  hydrocarbon  pyrolysis  and  oxidation  mechanism 
is  considerably  larger  than  the  previous  C4  model.  The 
C4  mechanism  included  447  reactions  among  127species, 
whereas  the  current  mechanism  contains  2498  reactions 
among  291  species. 

Fig.  3  illustrates  some  of  the  critical  steps  along  the 
pathway  from  a  parent  fuels  species  (e.g.  ethane)  to  a  de¬ 
posit  precursor  (e.g.  cyclopentene).  A  series  of  abstraction, 
beta- scission  and  addition/recombination  reactions  result 
in  a  stable  five-carbon  ring.  Under  the  appropriate  condi¬ 
tions,  similar  reaction  steps  can  continue  to  produce  larger 
molecules. 


2.3.  Solution  method 

The  governing  equations  may  be  a  system  of  ordinary 
differential  equations  (ODE)  or  differential- algebraic  equa¬ 
tions  (DAE),  depending  on  the  details  of  the  surface  chem¬ 
istry.  For  the  results  discussed  in  this  paper,  the  governing 
equations  are  a  set  of  ordinary  differential  equations.  In 
more  complex  surface-reaction  mechanisms,  which  incor¬ 
porate  surface  site  fractions  as  dependent  variables,  the 
equations  become  differential- algebraic  equations.  In  ei¬ 
ther  case,  the  DASSL  [18,19]  software  is  used  to  solve  the 
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Fig.  3.  Primary  reaction  pathway  from  ethane  to  cyclopentene. 


system.  DASSL  implements  a  variable  order,  variable  step 
size,  implicit  method  based  on  the  backward  differentiation 
formulas.  It  automatically  adjusts  the  step  size  to  control  lo¬ 
cal  truncation  error  according  to  a  user- specified  tolerance. 

In  addition  to  the  species  distributions  predicted  by  the 
kinetics  simulation,  the  local  equilibrium  state  is  also  deter¬ 
mined.  The  equilibrium  prediction,  which  constrains  pres¬ 
sure  and  temperature,  begins  with  the  kinetically  predicted 
species  distribution  at  a  particular  point  in  the  channel.  The 
free-energy  minimization  procedure  is  accomplished  with 
the  Equil  [20]  software,  which  uses  a  Chemkin-III  interface 
to  STANJAN-III  [21].  The  equilibrium  evaluation  considers 
all  the  species  that  are  included  in  the  kinetics  reaction  mech¬ 
anism.  Since  pure  carbon  (graphite)  is  not  included  in  the 
kinetics  it  is  also  not  included  in  the  equilibrium  prediction. 

2.4.  Parameters  and  boundary  conditions 

The  results  in  this  paper  are  based  on  a  10  cm  long  rectan¬ 
gular  channel  that  is  0.5  mm  high  and  1.0  mm  wide.  Assum¬ 
ing  1  mm  spacing  between  channels,  the  electrochemically 
active  area  per  unit  length  of  channel  is  a'e  =  2  mm.  Inlet 
velocities  of  40,  60  and  150  cm/s  are  used. 

The  temperature  profile  along  the  channel  is  uniform,  with 
calculations  reported  at  800  and  900  °C.  While  the  model  is 
capable  of  various  temperature  profiles,  these  effects  are  not 
explored  in  this  work.  The  feed  gas  is  95%  CH4/5%  C2H6, 
although  we  assume  that  only  methane  is  electrochemically 
oxidized.  Ethane,  commonly  found  in  commercial  natural 
gas,  participates  in  initiation  reactions  and  tends  to  accelerate 
molecular- weight-growth  chemistry. 

The  electrochemical  performance  of  a  particular  system 
requires  knowledge  of  the  system  architecture,  materials, 
charge-exchange  kinetics,  and  the  associated  overpotential 
losses.  Although  none  of  these  are  specifically  considered  in 
this  paper,  knowledge  of  the  current-density  profile  is  suffi¬ 
cient  to  determine  the  product  fluxes  into  the  anode  channel. 
We  consider  two  limiting  cases:  one  in  which  the  current 
density  is  uniform  along  the  channel  and  the  other  in  which 
the  current  density  has  first-order  dependence  on  the  fuel 
concentration.  The  order- dependence  parameter  y  in  Eq.  (2) 
is  0  for  a  uniform  current  density  and  1  for  a  current  den¬ 
sity  that  is  proportional  to  the  local  CH4  concentration.  The 
value  of  /g  is  chosen  to  provide  a  uniform  current  density  of 
600  mA/cm2  for  y  =  0  and  a  current  density  of  600  mA/cm2 
at  the  inlet  for  y  =  1.  These  limiting  electrochemical  as¬ 
sumptions  are  sufficient  to  achieve  our  objective  of  predict¬ 
ing  the  role  of  homogeneous  kinetics  in  coking  propensity. 


The  assumptions  are  clearly  not  sufficient  to  predict  cell  per¬ 
formance  generally. 


3.  Results 

By  varying  the  inlet  velocity,  and  hence  the  fuel  flow  rate, 
both  residence  times  and  composition  profiles  are  changed 
significantly.  Increasing  the  flow  rate  increases  the  flux  of 
fuel  down  the  channel  relative  to  the  flux  of  electrochemi¬ 
cal  reactants  and  products  through  the  MEA  active  area.  As 
shown  in  Fig.  4  for  y  =  0,  higher  fuel  flow  rates  lead  to  gen¬ 
erally  higher  CH4  concentrations  and  lower  H2O  and  CO2 
concentrations.  Residence  times  for  these  conditions  are  158, 
115  and  56  ms  for  inlet  velocities  of  40,  60  and  150  cm/s, 
respectively.  Though  some  CH4  reacts  homogeneously,  de¬ 
pletion  of  CH4  is  dominated  by  electrochemistry.  H2  and 
CO  mole  fractions  are  less  than  0.01  for  all  sets  of  param¬ 
eters  and  boundary  conditions,  supporting  the  assumption 
that  very  little  CH4  is  converted  homogeneously  to  these 
electrochemically  active  species.  Rather,  gas-phase  CH4  de¬ 
pletion  occurs  mainly  along  pyrolysis  reaction  pathways. 

For  a  cell  operating  at  900  °C,  several  cyclic  hydrocar¬ 
bon  species  are  produced  in  appreciable  quantities,  with 
cyclopentene  appearing  in  the  highest  concentration.  Fig.  5 
shows  the  mole  fraction  of  cyclopentene  for  the  various 
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Fig.  4.  Profiles  of  CH4,  H2O  and  CO2  mole  fractions  with  various  inlet 
velocities  at  900  °C  and  y  =  0. 
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Fig.  5.  Profiles  of  cyclopentene  mole  fraction  with  various  inlet  velocities 
at  900  °C  and  y  =  0. 


inlet  velocities  and  y  =  0.  The  lower  flow  rates  lead  to 
much  higher  levels  of  cyclopentene,  despite  the  elevated 
concentration  of  H2O  and  CO2.  This  is  mainly  attributed 
to  a  significant  increase  in  residence  time  and  correspond¬ 
ingly  an  increase  in  the  homogeneous  pyrolysis.  The  kinetic 
predictions  for  this  scenario  indicate  that  deposit  precursor 
formation  can  be  inhibited  by  reducing  the  residence  time, 
but  this  implies  a  decrease  in  fuel  utilization. 

Directly  comparing  the  kinetic  prediction  to  the  local  equi¬ 
librium  solution  provides  better  basis  for  understanding  the 
progression  of  molecular-weight  growth  chemistry  than  ki¬ 
netics  alone.  All  species  with  five  or  more  carbon  atoms  and 
cyclic  or  polycyclic  structure  are  considered  to  be  indicators 
of  a  propensity  to  form  deposits.  While  the  coking  precur¬ 
sors  produced  in  the  kinetic  model  are  typically  single-ring 
hydrocarbons,  the  equilibrium  model  tends  to  predict  high 
concentrations  of  multi-ring  aromatics.  To  compare  the  re¬ 
sults  of  the  two  models,  the  analysis  is  simplified  by  looking 
at  specific  groups  of  molecular  species.  The  fractional  car¬ 
bon  content  Fq  is  defined  as  the  ratio  of  the  number  of  car¬ 
bon  atoms  in  deposit  precursor  species  to  the  total  number 
of  carbon  atoms. 

As  seen  in  Fig.  6,  the  kinetic  and  equilibrium  predictions 
for  coking  propensity  are  substantially  different.  The  kinetic 
model  indicates  that  the  greatest  likelihood  for  deposit  for¬ 
mation  is  near  the  end  of  the  channel,  where  the  fuel  stream 
has  experienced  the  longest  residence  time.  The  fraction  of 
carbon  in  precursors  is  much  higher  for  cases  with  longer 


Fig.  6.  Profiles  of  fractional  carbon  content  in  deposit  precursors  for  both 
kinetic  and  equilibrium  models  with  various  inlet  velocities  at  900  °C  and 
y  —  0. 


residence  times.  At  higher  flow  rates,  the  equilibrium  pre¬ 
dictions  indicate  a  much  higher  propensity  to  form  deposits 
than  the  kinetic  model,  with  the  highest  fraction  of  carbon 
atoms  in  deposit  precursors  near  the  channel  inlet  where 
there  is  very  little  oxygen  in  the  system.  However,  note 
the  dramatic  drop  in  equilibrium  deposit  precursors  further 
downstream.  This  is  due  to  the  increased  oxygen  atom  con¬ 
tent  in  the  stream  as  more  H2O  and  CO2  enter  the  channel. 
Interestingly,  the  kinetic  model  and  equilibrium  calcula¬ 
tions  crossover  due  to  slow  rates  for  processes  involving 
molecular- weight  growth. 

Hydrogen  abstraction  from  water  by  radical  species  R, 
such  as 

H20  +  R  ^  RH  +  OH,  (R3) 

produce  OH  radicals  which  are  capable  of  oxidizing  large 
molecules.  However,  the  rate  of  (R3)  is  very  low  under  these 
conditions,  making  this  pathway  for  deposit  precursor  de¬ 
struction  ineffective.  Thus,  the  deposit  precursors  continue  to 
form,  but  at  a  slower  rate  than  near  the  inlet.  As  the  flow  rate 
increases,  the  shorter  residence  time  and  change  in  stream 
composition  cause  the  crossover  point  to  shift  downstream. 
The  homogeneous  kinetics  suggest  increases  in  deposit  pre¬ 
cursor  concentration  further  down  the  channel  in  spite  of  the 
increased  concentration  of  H2O  and  CO2.  One  might  be  able 
to  take  advantage  of  the  increased  oxygen  content  near  the 
outlet  by  introducing  an  active  catalyst  to  bring  the  system 
to  equilibrium  more  quickly. 

As  expected,  when  the  assumed  value  of  y  increases  from 
0  to  1,  the  rate  at  which  electrochemical  reactions  convert 
CH4  to  H2O  and  CO2  decrease  in  downstream  portions  of 
the  channel  as  CH4  is  depleted.  The  difference  in  electro¬ 
chemical  oxidation  rate  manifests  itself  in  the  relatively  high 
CH4  mole  fraction  and  relatively  low  H2O  and  CO2  mole 
fractions  shown  in  Fig.  7.  For  inlet  velocities  of  40,  60  and 
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Fig.  7.  Profiles  of  CH4,  H2O  and  CO2  mole  fractions  with  various  inlet 
velocities  at  900  °C  and  y  =  1. 
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Fig.  8.  Profiles  of  cyclopentene  mole  fraction  with  various  inlet  velocities 
at  900  °C  and  y  =  1. 

150  cm/s,  the  residence  times  are  increased  to  182,  124  and 
57  ms,  respectively.  The  increased  residence  time  and  de¬ 
crease  in  dilution  by  H2O  and  CO2  lead  to  significantly 
higher  levels  of  cyclopentene  in  the  fuel  stream.  Compar¬ 
ing  Figs.  5  and  8  shows  that  the  cyclopentene  mole  fraction 
more  than  doubles  when  the  functional  dependence  of  cur¬ 
rent  density  changes  from  zero  to  first  order. 

Fig.  9  shows  comparisons  between  kinetic  and  equilib¬ 
rium  solutions  when  the  electrochemistry  is  first  order  in  the 
fuel  concentration,  y  =  1.  Not  only  are  both  the  kinetic  and 
equilibrium  predictions  for  coking  propensity  higher  than 
that  for  y  =  0,  the  higher  order  dependence  leads  to  a  de¬ 
crease  in  the  current  density  as  CH4  is  depleted.  The  de¬ 
creased  current  density  causes  a  substantially  lower  average 
power  output. 

The  kinetic  model  predicts  a  very  strong  temperature  de¬ 
pendence,  as  seen  in  a  comparison  between  the  900  °C  calcu¬ 
lation  in  Fig.  6  and  an  analogous  calculation  at  800  °C  shown 
in  Fig.  10.  The  100  °C  temperature  reduction  decreases  the 
kinetically  predicted  Fq  by  several  orders  of  magnitude. 
The  equilibrium  calculation  changes  as  well,  but  with  a  less 
dramatic  decrease  in  coking  propensity.  While  decreasing 
the  operating  temperature  is  desirable  for  minimizing  de¬ 
posits  in  the  anode,  this  is  not  necessarily  a  practical  solution 
because  thermally  activated  electrochemical  processes  are 
less  effective  at  lower  temperatures.  The  optimal  operating 
temperature  must  consider  the  competition  between  coking 
propensity  and  electrochemical  performance. 


Fig.  9.  Profiles  of  fractional  carbon  content  in  deposit  precursors  for  both 
kinetic  and  equilibrium  models  with  various  inlet  velocities  at  900  °C  and 

y  =  f 
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Fig.  10.  Profiles  of  fractional  carbon  content  in  deposit  precursors  for  both 
kinetic  and  equilibrium  models  with  various  inlet  velocities  at  800  °C  and 
y  =  0. 

4.  Summary 

We  have  developed  a  modified  plug-flow  model  to  ana¬ 
lyze  both  the  homogeneous  kinetic  and  equilibrium  chem¬ 
istry  in  an  anode  flow  channel  of  a  DECO-SOFC.  This 
model  accounts  for  the  effects  of  electrochemical  oxidation 
of  hydrocarbon  fuel  on  both  transport  processes  and  reac¬ 
tion  kinetics.  While  the  reaction  processes  in  these  systems 
also  involve  complex  heterogeneous  chemistry,  focusing  on 
the  gas-phase  kinetics  and  equilibrium  provides  an  essential 
first  step  toward  understanding  fuel  chemistry.  While  the 
results  presented  here  are  for  dry  natural  gas,  the  model 
can  be  readily  extended  to  other  hydrocarbon  fuels  such  as 
propane  and  butane. 

The  homogeneous  kinetics  and  equilibrium  models  show 
an  interesting  crossover  in  coking  propensity.  Near  the  chan¬ 
nel  entrance,  where  there  is  relatively  little  oxygen  content, 
the  equilibrium  prediction  shows  high  coking  propensity. 
However,  owing  to  relatively  slow  reaction  rates,  the  kinetic 
predictions  show  small  concentrations  of  coking  precursors 
early  in  the  channel.  In  downstream  regions  of  the  channel, 
the  situation  changes.  With  increasing  levels  of  oxygen  (pri¬ 
marily  in  the  form  of  H2O  and  CO2),  chemical  equilibrium 
favors  volatile  carbon  compounds  and  deposit  formation. 
However,  again  owing  to  the  relatively  slow  homogeneous 
reaction  chemistry,  the  kinetics  prediction  continues  to 
show  increasing  levels  of  coking  precursors.  The  crossover 
point  depends  on  operating  temperature,  fuel  flow  rate  (res¬ 
idence  time),  and  the  rates  of  electrochemical  reaction.  The 
crossover  behavior  suggests  that  there  may  be  ways  to  use 
spatially  distributed  catalysts  that  selectively  promote  equi¬ 
librium  later  in  the  channels  where  equilibrium  suppresses 
coking.  Practical  implementation  of  this  approach,  how¬ 
ever,  would  likely  be  difficult.  The  functional  dependence 
of  current  density  on  fuel  concentration  strongly  influences 
the  tendency  to  foul  porous  anode  structures.  Higher  val¬ 
ues  of  y  will  increase  both  kinetic  and  equilibrium  coking 
propensity,  with  a  much  greater  increase  in  the  kinetic 
prediction.  The  higher  order  dependence  also  decreases 
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the  average  power  density  of  the  cell,  as  current  density 
decreases  with  fuel  depletion.  Temperature  also  has  a  sig¬ 
nificant  influence  on  molecular  weight  growth  chemistry, 
increasing  both  precursor  production  rates  and  the  equilib¬ 
rium  concentration  of  cyclic  species  and  multi-ring  aromatic 
hydrocarbons. 
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